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ABSTRACT 

The two-component, core-crust, model of a neutron star with homogenous in- 
ternal and dipolar external magnetic field is studied responding to quake-induced 
perturbation by substantially nodeless differentially rotational Alfven oscillations 
of the perfectly conducting crustal matter about axis of fossil magnetic field 
frozen in the immobile core. The energy variational method of the magneto- 
solid-mechanical theory of a viscoelastic perfectly conducting medium pervaded 
by magnetic field is utilized to compute the frequency and lifetime of nodeless tor- 
sional vibrations of crustal solid-state plasma about the dipole magnetic-moment 
axis of the star. It is found that obtained two-parametric spectral formula for the 
frequency of this toroidal Alfven mode provides fairly accurate account of rapid 
oscillations of the X-ray fiux during the fiare of SGR 1806-20 and SGR 1900+14, 
supporting the investigated conjecture that these quasi-periodic oscillations owe 
its origin to axisymmetric torsional oscillations predominately driven by Lorentz 
force of magnetic field stresses in the finite-depth crustal region of the above 
magnet ars. 

Key words: stars: individual (SGR 1900+14, SGR 1806-20); stars: neutron; stars: 
oscillations 



1. Introduction 

The quasi-periodic oscillations (QPOs) recently detected in the X-ray fiux of fiaring 
SGR 1806-20 and SGR 1900+14 and interpreted as a manifest of quake-induced torsional 
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vibrations of underlying magnetar (Israel et al 2005, Watts & Strohmayer 2006) are among 
the topical themes in current development of the astroseismology of this subclass of neutron 
stars (e.g. Glampedakis et al 2007; Samuelssen & Andersson 2007; Levin 2007, van Hoven & 
Levin 2008, Lee 2007, 2008; Sotani et al 2007, 2008; Vavoulidis et al 2008) whose bursting and 
quiescent electromagnetic activity is dominated, as is commonly believed, by frozen-in fossil 
magnetic field of extraordinary intensity (e.g.. Woods & Thompson 2006, Mereghetti 2008). 
Motivated by this interest, in (Bastrukov et al 2007, 2008, 2009) several possible scenarios 
of post-quake relaxation of magnetar has been investigated with emphasis on regime of 
nodeless torsional vibrations about the dipole magnetic-moment axis of these seismically 
active neutron stars. In particular, in (Bastrukov et al 2007b, 2008), a case of the elastic- 
force-driven nodeless shear oscillations, both torsional - otf and spheroidal - oSe, locked in 
the finite-depth crust AR has been studied in some details. It was found that obtained 
two-parametric spectral equations for the frequency of pure elastic spheroidal and torsional 
modes as a function of multipole degree i of nodeless shear oscillations of the form (Bastrukov 
et al 2008) 



provide proper account of the low-frequency QPOs frequencies, supporting £-po\e identifi- 
cation proposed in works (Watts & Strohmayer 2007, Samuelssen & Andersson 2007) and 
it was shown that the unidentified before lowest overtones of QPOs in light-curve of flar- 
ing SGR 1806-20 with frequencies 18 and 26 Hz can be interpreted as manifest of dipole, 
i — 1, spheroidal and torsion nodeless elastic shear vibrations, exhibiting features generic to 
Goldstone soft modes. 

On the other hand, one can perhaps cast doubt on such interpretation of QPOs as being 
produced by seismic vibrations driven by solely elastic restoring force, because it rests on a 
poorly justifiable assumption about dynamically passive role of ultra-strong magnetic field 
defining only direction of axis about which the stellar matter undergo torsional nodeless 
oscillations and remaining unaltered in the process of post-quake relaxation of magnetar. 
Instead, it seems more plausible to expect that explosive liberation of magnetic field energy 
resulting in the crust fracture by relieved magnetic field stresses - magnetic breaking associ- 
ated with disruption of the magnetic field lines - must be followed, after the magnetic field 
lines reconnection, by the Lorentz-force-driven Alfven differentially rotational oscillations of 
the perfectly conducting star matter about axis of frozen-in homogeneous magnetic field, 
either in the entire volume or in the peripheral finite-depth crust. This line of argument has 
been taken into account in the magneto-solid-mechanical model of global torsional nodeless 
vibrations about the dipole magnetic moment axis of neutron star driven by the joint action 



Uj{oSe) = a;s(£,CUe, A), Uj{oti) = UJt{C-,UJe, 



£>1, 



(1) 




(2) 
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of Hooke's force of elastic shear stresses and Lorentz force of magnetic field stresses (Bas- 
trukov et al 2009). The computed in this latter work the two-parametric frequency spectrum 
of the form 

u;{ote)^u;t{i,u;,,P), /3 = ^ = £ > 2, (3) 

= = ? = "^ = 4^ 

was found to be consistent with data on QPOs with frequencies u = 00/271 from the range 
30 < z/ < 200 Hz interpreted as being produced by these global seismic vibrations of multipole 
degree I from interval 2 < £ < 20, but unable to explain the above mentioned lowest overtones 
in the QPOs spectrum of SGR 1806-20. Also, from the viewpoint of both outlined models 
the nature of restoring force responsible for high-frequency QPOs with v — 626 and v — 1837 
Hz remains uncertain. 

In this paper, continuing investigations reported in (Bastrukov et al 2007, 2008, 2009), 
we focus on different, to the above mentioned scenarios, case of post-quake vibrational re- 
laxation of a neutron star. Specifically, we consider torsional nodeless Alfven oscillations of 
the crustal solid-state plasma, that is, seismic vibrations restored by solely Lorentz force of 
magnetic field stresses and excited in the peripheral finite-depth seismogenic layer of neutron 
star model with homogeneous internal and dipolar external magnetic field. In section 2, a 
brief outline is given of equations of magneto-solid-mechanics appropriate for the perfectly 
conducting viscoelastic continuous medium pervaded by magnetic field. In section 3, the 
spectral formulae for the frequency and lifetime of this toroidal Alfven vibration mode are 
obtained followed by discussion of inferences that can be made from considered scenario 
of seismic vibrations regarding electrodynamic properties of the neutron star matter with 
homogeneous configuration of the frozen-in magnetic field. The relevance of considered aster- 
oseismic model to the detected QPOs in the X-ray flux during the flare of SGR 1806-20 and 
SGR 1900-1-14 is assessed in section 4. The newly obtained results are briefly summarized 
in section 5. 



2. Governing equations of magneto-solid-dynamics of viscoelastic perfectly 
conducting stellcir matter in the presence of frozen-in magnetic field 

The equation of magneto-solid-dynamics describing mechanical distortions of the per- 
fectly conducting continuous medium of extremely high electric conductivity threaded by 
magnetic fleld are 

pili = Vfe Tik Vjfc TTjfe (5) 



-4- 



where Ui is the field of material displacements - the basic variable of the Newtonian elastic 
dynamics of solid material continuum, 

Tik = ]- [B:6Bk + Bk6Bi - B,6BM , (5B = V x [u x B] (6) 

4:71 

stands for the tensor magnetic field stresses (Franko et al 2000, Bastrukov et al 2009) and 

TTik = 2r]Uik, iiik = ^[VjMfc + VfcMj] (7) 

is the tensor of viscose stresses vTjfc linearly related to the rate of shear strains iiik with rj being 
the shear viscosity of stellar matter. This last constitutive equation expresses Newtonian law 
of viscosity which holds for both liquid and solid aggregate states of continuous medium. The 
equation of energy balance in the process of compression-free Alfven vibrations highlighting 
its shear character reads 

d f fru? 



J ^ dV = - J Tik Uik ^ J '^ik ^i'' '^i^k = ^kUk = 0. (8) 

This equation is obtained by scalar multiplication of equation of magneto-solid-mechanics 
([5]) with Ui and integration over the star volume. To compute the frequency and lifetime 
of Alfven oscillations, we take advantage of Rayleigh's energy variational method which 
assumes the following separable representation of fluctuating field of material displacements 

Ui{r,t) = ai{r)a{t) (9) 

with aj(r) being the time- independent field of instantaneous displacement and amplitude a{t) 
describes temporal evolution of fluctuations. The tensors of fluctuating Maxwell's magnetic 
stresses and Newtonian viscous stresses are too written in similar separable form 

Tik{r,t) = [fik{r) - ^fjj{r)Sik]a{t), (10) 

f,,(r) = -^[5,(r) bkir) + 5fe(r), 6,(r)], (11) 
47r 

k{v) = Vfc [a,(r) Bkir) - afe(r) 5,(r)], (12) 

Tiikir) = 2r]aik{r)a{t), aifc(r) = -[Viafc(r) + Vfcai(r)]. (13) 

The idea of above separable representation of fluctuating variables is that on substituting 
(!9|)- (fT3|) in (jH]), this latter equation is reduced to equation for a(t) having the well-famihar 
form of equation of damped oscillations 

Ma + 'Dce + ICa = 0, a(t) = ao exp(-t/r) cos(rit), (14) 
n' = u^[l~{urr], u^ = ^, r=^. (15) 
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with the integral parameters of the inertia M., viscous friction T> and stiffness /C given by 

J p{r)ai{r)ai{r)dV, (16) 

= ^ y [Bi{r)(3k{r) + 5fc(r)A(r)] [VM^) + V^a^lr)] dV, (17) 

^Wt'^^"*^ W + [Viafe(r) + Vfeai(r)] dV. (18) 

The above expounded equations are general in the sense that they can be used for computing 
frequency uo of both the even-parity poloidal and odd-parity toroidal Alfven vibrational 
modes (Bastrukov et al 1999) and the time of viscous damping of nodeless oscillations in 
a neutron star models in which the density, the magnetic field and the shear viscosity are 
considered to be arbitrary functions of position. 

As was noted by McDcrmott ct al (1988), the characteristic feature of solid- mechanical 
treatment of the neutron star astcroscismology, having many features in common with the 
one adopted in geoseismology (e.g. Lapwood & Usami 1981, Lay & Wallace 1995, Aki & 
Richards 2003), is that it deals with quake-induced shear vibrations about hydrostatic equi- 
librium which is presumed to exist, but with no reference to explicit form of equation of 
gravitationally equilibrium state which is normally considered as the free from anisotropic 
stresses state. The anisotropic stresses of different physical nature (elastic, magnetic or grav- 
itational) defining restoring force of shear vibrations are regarded to come into play only after 
quake-induced perturbation. By thus relinquishing the condition of hydrostatic equilibrium 
of magnetic solid star, we gain enormous freedom in the choice of possible configurations of 
magnetic field inside and outside the star with the only constraint that this equilibrium field 
must obey the equations of magneto-statics. The remainder of this paper, which is sequel 
to (Bastrukov et al 2007b, 2009), considers the neutron star model with homogeneous inside 
and dipolar outside magnetic field for which this condition is of course met. 



3. Toroidal Altven mode of nodeless differentially rotational oscillations in the 
peripheral finite-depth seismogenic crust of a neutron star 

The application of the magnetic solid star model under consideration to neutron stars 
rests on the common today behef that the frozen in the stellar material magnetic field is 
fossil (e.g. Flowers & Ruderman 1977, Bhattacharya & van den Heuvel 1991, Chanmugam 
1992, Spruit 2008). Contrary to the main sequence (MS) stars whose magnetic fields are 
generated by permanent fiows of fiuid in the self-exciting dynamo processes, in the neutron 
stars there are no nuclear energy sources to support convection, so that the fossil magnetic 
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field is considered to be frozen in tlie super dense solid material continuum which before 
starquake remains immobile being highly compressed by immense pressure of self-gravity 
which is counterbalanced by the degeneracy pressure of Fermi-gas of relativistic electrons in 
the star crust and non-relativistic neutrons in the core. The magnetic field pressure con- 
structively contributes to the total internal pressure, however its absolute value is much less 
than the degeneracy pressure, so that magnetic pressure plays insufficient role in hydrostatic 
equilibrium of pulsars and magnetar^, but may play substantial, if not decisive, part in the 
process of their post-quake relaxation (Bastrukov et al 2009). 

In spherical polar coordinates, the components of internal homogeneous field directed 
along the polar axis are given by 

Br = B cos 9, Be = -B sin 9, B^ = (19) 

and external dipolar configuration is described by the field B = V x A, where A = [0, 0, A^j) = 
m^/r^] is the vector potential with the standard parametrization of the dipole magnetic 
moment m^ = {1/2)BR^ of the star of radius R and by B is understood the magnetic field 
intensity at the magnetic poles of the star, B = Bp, as is illustrated in Fig.l. 

The material displacements in the crust undergoing non-radial differentially rotational 
oscillations about magnetic axis of the star are described by the nodeless toroidal field of the 
form (Bastrukov et al 2007b, 2008b) 

u(r, t) = at(r) a(t), at(r) = [0(r) x r] (20) 
0(r) = Vx(r), V\{r) = (21) 
xir) = [Aer' + Ber''-']Pe{C), C = cos9 (22) 

The arbitrary constants Ae and Be are eliminated from the following boundary conditions 

U^\r=R. = 0, U^\r=R = X R]^, $ = V.=rP,(C), Vr=R = (^0, -^^-^^^ . (23) 

The leftmost equation is the no-slip condition on the core-crust interface, r = R^, implying 
that the amplitude of differentially rotational oscillations is gradually decreases as depth of 



^At normal nuclear density p — nirin — 2.8 10^^ g cm"^ with m being the mass of neutron, the pressure 
of degenerate Fermi-gas of non-relativistic neutrons Pf = Kmru^^^ « 10^^ dyn cm~^. The typical for 
magnetars magnetic fields 10^* < B < 10^^ Gauss, so that the magnetic field pressure, Pb — B'^/{8tt), 
falls in the interval 10^^ < Pm < 10"^^, that is, at least two order of magnitude less than Pp. For this 
reason we do not take into account the Chandrasekhar- Fermi effect of magnetic flattening of spherical shape 
of star at its magnetic pole which has first been established for homogeneous magnetic field of dipolar 
external configuration (Chandrasekhar & Fermi 1953) and confirmed by Ferraro (1954) for a star model with 
nonhomogeneous poloidal internal and dipolar external field. 



Fig. 1. — The dipolar external and homogenous internal magnetic field of the two- 
component, core-crust, model of the neutron star under consideration. 



seismogenic layer is increased and turns into zero on the surface of core. The condition on 
the star surface, r = R, is dictated by symmetry of the general toroidal solution of the vector 
Laplace equation as well as the fact that resultant expression of mass parameter of torsional 
oscillations must lead to correct expression for the moment of inertia of a rigidly rotating star 
(Bastrukov et al 2008). Explicitly, the above boundary conditions are described by coupled 
algebraic equations 



AiR'-' + BeR^ 
AiR^ + B,R-^-^ 



-1-2 



0, 



R 



whose solutions read 

Ai = Mi, Be 



-Me Rf^^ , Afe = -^eTTZT^e+T ■ 



(24) 
(25) 



(26) 



The spherical components of the nodeless toroidal field of instantaneous displacements are 
given by 

Be 



0, ae = 0, 



Aer' + 



'I _ ^2n1/2 '^-P^(C) 



(27) 



The snapshot of quadrupole, £ = 2, and octupole £ = 3, overtones of these axisymmetric 
differentially rotational vibrations about the dipole magnetic-moment axis of a neutron star 
is pictured in Fig.2. 
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Fig. 2. — Quadrupole and octupole overtones of the nodeless differentially rotational, tor- 
sional, vibrations of crustal matter about the dipole magnetic moment axis of neutron star. 



The mass parameter of torsional nodeless vibrations, computed in the approximation of 
uniform density of crustal matter, reads 



M 

mt{i, A) 
A 



4:71 pR^ 



+ 1) 



(2£+l)(2£ + 3) 



mt{i, A), 



(2^ 



(1 _ a2^+1)-2 

Rr 



(2£ + 3)A 



2e+i 



+ 



i2£ + l] 
2£-l 



-A 



2e+3 



2i + 3 
2i-l 



^2(2£+l) 



R 



< A < 1. 

The integral parameter of viscous friction D is given by 

d,(£,A) = (l-A2^+i)-i 



V{ote)=A7rr]R'^-^^-^dt, 



21+ I 



1 - 



{^ + 2) .2^+1 



(29) 



(30) 



For the lifetime of torsional oscillations damped by shear viscosity r/ of crustal material we 
obtain 



2t„ 



mtji, A) 

(2£ + 3)(£-l) dt{e,X) 



R^ 



V 



V 

V = —. 

P 



(31) 



In Fig. 3, the fractional lifetime is plotted as a function of multipole degree C. of torsional 
nodeless oscillations for indicated values of fractional depths h of crustal seismically active 
region of the neutron star. 

At this point it may be appropriate to note that replacing, in the last equation for 
T>[ot£), the coefficient of shear viscosity rj by shear modulus /i of crustal matter, the resultant 
equation becomes identical to the expression for stiffness Ke{ote) of the elastic- force-driven 
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Fig. 3. — The fractional lifetime of torsion nodeless oscillations in the neutron star crust as 
a functions of multipole degree I and indicated values of the fractional depth h of peripheral 
seismogeneic zone. 

nodeless oscillations (Bastrukov et all 2007b, 2008) whose general integral form is given by 

^e = \j Ai(r)[Viafc(r) + Vfeai(r)] [Viafc(r) + Vfeai(r)] dV (32) 

where /i(r) stands for the shear modulus. Thus, in the standard two-component, core- 
crust, model of quaking neutron star (Franco et al 2000) presuming homogeneous density 
profile for crustal matter p and constant values of transport coefficients of shear viscosity 
rj and shear elasticity /i one has Ke{ote) = {fi/r])V{Qt£); the frequency of elastic nodeless 
vibrations is estimated by ratio u'^ = ICe/Ai (Bastrukov et all 2007b, 2008). In the limit, 
A = (Rc/R) 0, the analytic expression for the damping time, r = 2A4/V, of torsional 
nodeless oscillations in the crust is reduced to that for hfetime of global torsional nodeless 
oscillations in the entire volume of homogeneous neutron star model (Bastrukov et al 2003, 
2007a) 

^^'''^ = (2£ + 3)(£-l)- 

The practical significance of outlined equations is that they can be used for the study of 
nodeless torsional vibrations of more wide class of celestial objects having solid seismogenic 
mantle of finite depth, such as Earth-like planets (Lay & Wallace 1995, Aki & Richards 2003) 
and white dwarf stars whose asteroseismology is the subject of intense current research (e.g. 
Winget & Kepler 2008). 

Computation of integral for the spring constant /Cm of the magnetic Lorentz restoring 
force of non-radial Alfven nodeless differentially rotational vibrations in the crust yields 
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(l_A2^+i)-2Jl + 



(£2-l)(2£ + 3) 



As a result we arrive at the spectral equation for the frequency as a function of multipole 
degree of nodeless Alfven oscillations 



{f - 1) 



va 



2i + 3 
2i-l 
B 



Va 



m,(£,A)' 



(35) 
(36) 



which depends, as one sees, upon two parameters - Alfven frequency, va — {1/2t:)va/ R 
and the fractional depth, h — AR/R = 1 — A, of seismogenic layer. In Fig.4 and Fig.5, 
the fractional frequencies and periods of toroidal Alfven vibration mode as a functions of 
multipole degree i with indicated values of fractional depth of the seismogenic layer h are 
plotted, showing that the higher i the larger is the frequency and the shorter is the period 
of Alfven nodeless vibrations. 

One sees that the lowest overtone of torsion Alfven oscillations entrapped in the crust 
is of dipole degree, i = 1. However, in the limit A = {Rc/ R) — > [or, in other terms, 
when h ^ 1 and kl^{i , X) / mt{i , \) — * 1], the case when entire volume of the star sets in 
the oscillations, the above two-parametric spectral formula is reduced to the one-parametric 
frequency spectrum of global Alfven torsional vibrations (Bastrukov et al 2009) 



z/(o4) 



l^A 



{f - 1) 



21 + 3 
21-1 



.1/2 



l^A 



UJa_ 

27r' 



uja 



RB^ 
~3M 



M 



pR^ (37) 



with the lowest overtone of quadrupole degree, 1 — 2. This suggests that dipole vibration 
is Goldstone's soft vibrational mode whose most conspicuous property is that the mode 
disappears (the frequency tends to zero) when one of parameters, in our case A, of vibrating 
system turns into zero (Bastrukov et al 2008). 

Physically, the considered solid-mechanical model can be invoked to analysis of quasi- 
periodic oscillations of the electromagnetic emission from neutron star with approximately 
homogenous internal and dipolar external magnetic field when quake-induced perturbation is 
fairly week, that is, does not penetrate into the deep interior of the star setting in oscillations 
only its peripheral solid mantle of finite depth. Also, the model is appropriate for a case 
when the electro dynamic properties of the core material are incompatible with conditions of 
existence of Alfven oscillations, one of the major of which is the extremely large (effectively 
infinite) electrical conductivity of stellar material (e.g., Chandrasekhar 1961, Mestel 1999). 
For the solid crustal matter, which is normally considered as a metal-like solid-state plasma 
composed of nuclei embedded in the degenerate Fermi-gas of relativistic electrons (whose 
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Fig. 4. — Fractional frequency and period of nodeless torsional magneto-solid-mechanical 
oscillations, toroidal Alfven mode - qO^, entrapped in the neutron star crust as functions 
of multipole degree I at different values of the fractional depth h of peripheral seismogenic 
crust. The value h = 1 corresponds to global torsional oscillations excited in the entire 
volume of the star. Here ua = i^A/27r, where ua = va/R with va = B / y/Anp being the 
velocity of Alfven wave in crustal matter of density p and Pa = Iti /uja- 



pressure prevents the crust form gravitational contraction) this condition is certainly fulfilled 
(e.g. Chanmugham 1992). However, this is may not be the case for much denser matter of 
the neutron star cores, because its micro-composition is dominated, according to canonical 
notion of the neutron star, by degenerate Fermi-matter of non-relativistic neutrons (whose 
pressure is responsible for the neutron star gravitational equilibrium) - the electrically neutral 
Fermi-particles endowed with intrinsic spin magnetic moments. This latter property of 
neutrons suggests that the degenerate Fermi-material of the neutron star cores can be brought 
to gravitational equilibrium in the magnetically ordered state of paramagnetic permanent 
magnetization caused by Pauli's mechanism of field-induced polarization of spin magnetic 
moment of neutrons along the axis of frozen in the core fossil magnetic field. It is just the 
case when the non-conducting core matter is capable of sustaining homogeneous magnetic 
field, but is unable to support Alfven oscillations. The structure of such paramagnetic 
neutron star model, which is extensively discussed in (Bastrukov et al 2002a, 2002b, 2003) 
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Fig. 5. — Fractional frequency of nodeless torsional Alfven oscillations of indicated overtone 
£ as a function of the fractional depth h of peripheral seismogenic layer. The vanishing 
of dipole overtone in the limit oi h —>■ 1, the case when entire mass of neutron star sets 
in torsional oscillations about dipole magnetic moment of the star value exhibits behavior 
which is typical for Goldstone's soft modes. 



is schematically pictured in Fig. 6. 

It worth noting that internal constitution of such a star does not exclude the existence 
of a fairly thin inner crust which is predominately composed of boson matter whose micro- 
constituents are paired neutrons with average density several times less than the normal 
nuclear density, 0.1 — 0.5pN- Such a matter permits Bose-condensate state having some 
features in common with that generic to weakly compressed superfluid liquids. Understand- 
ably that such boson fraction of neutron matter, however, can only insufficiently contribute 
to the total mass budget of neutron star, because Bose-matter is unable to withstand the 
pressure of self-gravity, contrary to the degenerate neutron Fermi-matter whose pressure of 
degeneracy counterbalancing the self-gravity pressure is central to the very notion of neutron 
star. The physical importance of thin inner crust is, most likely, in that it operates like a lu- 
bricant facilitating differentially rotational shear of conducting crust relative to much denser 
paramagnetic core. From the view point of this model, the coupling between core and crust 
is one and the same of physical nature of magnetic cohesion as that operating between sam- 
ples of terrestrial metal and permanent magnet which is mediated by the lines of magnetic 
force. The idea of decisive role of magnetic core-crust cohesion seems to be in line with 
the starquake model of the magnetar burst (Woods and Thompson 2006) which is though 
of as disruption of magnetic field lines in the process of explosive release of magnetic field 
stresses fracturing the neutron star crust. The presented line of argument (while admittedly 
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PARAMAGNETIC NEUTRON STAR 




Fig. 6. — The internal constitution of two-component, core-crust, model of paramagnetic 
neutron star model. The core is considered as a poorly conducting permanent magnet com- 
posed of degenerate Fermi-gas of non-relativistic neutrons in the state of paramagnetic satu- 
ration caused by Pauli's paramagnetism - field-induced alignment of spin magnetic moment 
of neutrons along the axis of uniform internal and dipolar external magnetic field frozen in 
the star on the stage of gravitational collapse of its massive progenitor. A highly conducting 
metal-like material of crust, composed of nuclei suspended in the degenerate Fermi-gas of rel- 
ativistic electrons, is regarded as magneto-active electron-nuclear solid-state plasma capable 
of sustaining Alfven oscillations. 

speculative, but no more than any other considerations touching upon the internal structure 
of neutron stars) shows what inference can be made regarding electromagnetic properties of 
super dense stellar matter from asteroseismic model of Lorentz-force-driven torsional Alfven 
nodeless oscillations entrapped in the finite-depth crust, provided of course that the above 
derived two-parametric spectral formula properly describe observational data on oscillating 
electromagnetic emission from quaking neutron star. In the next section we assess this con- 
jecture by applying the derived theoretical spectrum to modal £-po\e classification of QPOs 
in the X-ray flux during the flare of SGR 1806-20 and SGR 1900+14. 
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4. QPOs in the X-ray flux during the flare of SGR 1806-20 and SGR 1900+14 
from the viewpoint of above asteroseismic model 

The physics behind interconnection between oscillating lines of magnetic field frozen in 
the star and quasi-periodic oscillations of detected electromagnetic fiux has been recognized 
long ago (e.g. van Horn 1980, Blaes et al 1988) - the quake-induced perturbation excites 
coupled oscillations of perfectly conducting matter and the lines of frozen-in magnetic field 
which outside the star operates as transmitters of beam of charged particles producing coher- 
ent radiation from the neutron star magnetosphere, so that the frequency of Alfven seismic 
oscillations is considered to be equal to the frequency of quasi-periodic variations of the beam 
direction; this asteriseismic argument regarding physical nature of quasi-periodic oscillations 
is applicable to both radio emission from glitching pulsars and high-energy emission from 
bursting magnetars. 

240 
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Fig. 7. — Theoretical fit of the QPOs frequency in the X-ray fiux during the fiare of SGRs 
1806-20 and 1900+14 by the obtained two-parametric spectral equation for the frequency of 
nodeless torsional magneto-solid-mechanical oscillations entrapped in the peripheral finite- 
depth crustal region of a neutron star with homogeneous internal and dipolar external mag- 
netic field. 

As was noted in our previous study of global Alfven oscillations in the entire volume of 
neutron star with homogeneous internal and dipolar external field (Bastrukov et al 2009), 
it is impossible reconcile one-parametric frequency spectrum i^(oo^) plotted as a function of 
multipole £, given by equation (!37|) . with the observed QPOs frequencies (measured in Hz), 
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which, according to compilation by Sotani et al (2008), for SGR 1806-20 are 18, 26, 29, 92, 
150, 625, 1837, and for SGR 1900+14 these are 28, 54, 84, 155. Bearing this in mind, one 
can attempt to fit the detected QPOs on the basis of derived here two-parametric spectral 
equation for torsional Alfven nodeless oscillations entrapped in the finite-depth peripheral 
crust. The adjustable parameters of such a fit are the Alfven frequency i^a and the fractional 
depth of seismogenic layer h. The suggested i-pole identification is summarized in Figs. 7 
and 8. One sees that at indicated in this figures parameters, the model adequately regain the 
general tendency in the data on QPOs frequencies when i are specified by integer numbers 
clearly recognizable in ticks of x-axis of Fig. 7. Deserved for special comment is that the 
low-frequency QPOs in data for SGR 1806-20, are fairly reasonably interpreted as dipole and 
quadrupole overtones of torsional magnetic Alfven nodeless vibrations, namely, z/(oai) = 18 
and ^(00-2) = 26 Hz. The high-frequency kilohertz vibrations with 627 Hz and 1870 Hz in 
the QPOs spectrum of SGR 1806-20 highlighted in Fig.8 are unambiguously specified as 
high-multipole overtones, namely, as z^(oCf^=42) = 627 Hz and ^(0*^^=122) ~ ^ 

870 Hz. 

5000 M I 




Fig. 8. — The same as Fig. 7, but for only SGRs 1806-20 with account for high-frequency 
QPOs. 

From the indicated in these figures values of z/^ = {1/2tt)va/R it follows that the 
average speed of Alfven wave va in the crustal matter falls in the range 10^ < va < 10*^ cm 
sec~^ and, thus, less than the average speed of elastic shear wave whose generally accepted 
value is q ~ 10^ cm sec~^. Also noteworthy is that the existence of free vibrations with 
frequency u = 2txv is determined by condition ujt » 1, as follows from equation (1151) . which 
is too fulfilled; but this of course does not mean that the conversion of energy of Alfven 
seismic crustal vibrations into heat (thermal emission) is the dominant mechanism of their 
exponential damping. All the above suggests that considered model of torsional nodeless 
Alfven oscillations locked in the crust, provides somewhat better account of the data on 
the QPOs frequencies, as compared to the model of global torsional nodeless vibrations of 
magnetar under combined action of elastic Hooke's and magnetic Lorentz forces (Bastrukov 
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et al 2009). 

5. Summciry 

The neutron stars, both pulsars and magnetars, are seismically active compact objects of 
the end point of the evolutionary stellar track. They are, as is widely beheved, to be formed 
in the magnetic-flux-conserving supernova collapse of the cores of massive stars coming 
into existence, as sources of pulsed electromagnetic emission, with a fairly complex internal 
constitution (e.g. Weber 1999, Bisnovatyi-Kogan 2002, Lattimer & Prakash 2007). The past 
two decades have seen increasing recognition that valuable tool of looking into the neutron 
star interior provides asteroseimology. The explosive interest to current development of this 
domain of pulsar astrophysics has been stimulated by the above mentioned discovery of 
quasi-periodic oscillations during the giant X-ray flares of SGR 1806-20 and SGR 1900-1-14, 
which has also prompted extensive investigations of radiative processes caused by ultra- 
strong magnetic fields of magnetars (e.g. Heyl & Hernquist 2005, Harding & Lai 2006, 
Zhang et al 2007, Timokhin et al 2008, Bo Ma et al 2008, Rea et al 2008) and processes of 
their cooling (Yakovlev et al 2008). 

Many agree today that detected QPOs can be explained on the basis of asteroseismic 
models of a magnetar undergoing quake-induced torsional oscillations about its the dipole 
magnetic-moment axis. Adhering to this attitude and taking into account that the main 
energy source of bursting high-energy emission of magnetars is the immense store of the 
magnetic field energy, we have set up and investigated one of such models interpreting the 
above QPOs in terms of axisymmetric Alfven nodeless oscillations of the crustal solid-state 
plasma about axis of homogeneous internal and dipolar external magnetic field frozen in the 
motionless core. As is demonstrated in Figs. 7 and 8, the obtained two-parametric spectral 
equation of this model provides proper account of general trends in data on QPOs frequencies 
as well as quite reasonable their £-po\e identification. The main argument in favor of such 
interpretation is fairly accurate description of both the low-frequency (18 and 26 Hz) and 
high-frequency (627 and 1837 Hz) QPOs in X-ray fiaring SGR 1806-20, which is presented 
in Fig.8, supporting the underlying this model conjecture about dominate role of Lorentz 
force of magnetic field stresses in maintaining the nodeless torsional oscillations of perfectly 
conducting solid-state plasma entrapped in the peripheral finite-depth seismogenic zone of 
magnetar. 

The physically meaningful framework for discussion of considered scenario of post-quake 
relaxation of magnetar provides the two-component, core-crust, model of quaking paramag- 
netic neutron star (Bastrukov et al 2002a, 2002b, 2003) which attributes the core matter 
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incapability of supporting Alfven differentially rotational shear oscillations to poor electrical 
conductivity of central region of the star, considering the magnetar core as a motionless 
spherical bar magnet composed of degenerate Fermi-matter of non-relativistic neutrons with 
spin magnetic moments polarized along the uniform fossil magnetic field frozen in the star 
on the stage of catastrophic contraction of its massive progenitor. Together with this, one 
cannot of course state that such an image of the quaking magnetar interior is unique, that is, 
that the above QPOs cannot be explained by alternative conjectures about electrodynamic 
properties of neutron star material, restoring forces and volumetric character of torsional 
oscillations (e.g. Bastrukov et al 2007b 2008, 2009). One more of such uninvestigated al- 
ternatives is the regime of nodeless global differentially rotational vibrations of neutron star 
with non-homogenous internal and dipolar external magnetic field presuming, thus, that 
stellar material in the total volume of magnetar possesses properties of perfect conductor. 
We shall return to such a model of quaking neutron star in a forthcoming paper. 

This work is a part of projects on investigation of variability of high-energy emission 
from compact sources supported by NSC of Taiwan, grant numbers NSC-96-2628-M-007- 
012-MY3 and NSC-97-2811-M-007-003. 
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